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Abstract 

We calculate the mean nucleation distance, dnuc^ in a first order cosmological 
quark-hadron phase transition. For homogeneous nucleation we find that 
self-consistent inclusion of curvature energy reduces dnuc to ^ 2cm. However, 
impurities can lead to heterogeneous nucleation with dnuc of several meters, 
a value high enough to change the outcome of Big Bang nucleosynthesis. 
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I. INTRODUCTION 



A number of authors have suggested that inhomogeneities in the baryon number 
density were produced in the cosmological quark-hadron phase transition and perhaps even 
persisted to the time of nucleosynthesis, which could alter the abundances of light elements. 
For this to happen, a few criteria have to be met, though. First of all, only a first or- 
der quark-hadron transition seems capable of generating large baryon number fluctuations. 
Whether the transition is first order is still an unsettled question (see e.g. Ref. 0). Sec- 
ondly, the mean distance between high and low baryon density regions, denoted /, has to 
be larger than approximately one meter comoving at T = lOOMeV [Q, if inhomogeneous 
nucleosynthesis results are to differ from the standard homogeneous nucleosynthesis results 
(otherwise, neutron and proton diffusion will erase the inhomogeneities prior to nucleosyn- 
thesis). Other important parameters are the average density contrast between high and 
low-density regions, the average volume fraction of high-density regions, and the baryon to 
photon ratio. 

We assume that the transition is first order and focus in this paper on determination of 
the most important of the parameters, the mean distance /. 

A first order quark-hadron transition in the expanding Universe can in general terms 
be described as follows. As the color deconfined quark-gluon plasma cools below the crit- 
ical temperature, ~ lOOMeV, it becomes energetically favorable to form color confined 
hadrons (primarily pions and a tiny amount of neutrons and protons due to the conserved 
net baryon number). However, the new phase does not show up immediately. As is char- 
acteristic for a first order phase transition, some supercooling is needed to overcome the 
energy expense of forming the surface of bubbles of the new hadron phase. When a hadron 
bubble is nucleated, latent heat is released, and a spherical shock wave expands into the 
surrounding supercooled quark-gluon plasma. This reheats the plasma to the critical tem- 
perature, preventing further nucleation in a region passed by one or more shock fronts 
(for the range of parameters used in the present investigation, bubble growth is described 
by deflagrations, where a shock front precedes the actual phase transition front. However, if 
the shock fronts are very weak the quark-gluon plasma may not be reheated enough to suffi- 
ciently prevent new bubble nucleations in such regions 0. This reduces the mean distance. 
We will neglect this possibility, because of the very small amounts of supercooling we deal 
with here). The nucleation stops when the whole Universe has reheated to Tc. This part of 
the phase transition passes very fast, in about 0.05/is, during which the cosmic expansion is 
totally negligible. After that, the hadron bubbles grow at the expense of the quark phase 
and eventually percolate or coalesce; almost in equilibrium at Tc. The transition ends when 
all quark-gluon matter has been converted to hadrons, neglecting possible quark nugget 
production 

How much the system has to supercool before the first hadron bubbles show up deter- 
mines the distance between nucleation sites. This distance is called the (comoving) mean 
nucleation distance, dnuc- It is also roughly the distance between the shrinking, high-baryon- 
density quark droplets. / is half that distance, and when we take into account that the 
Universe expanded about 50% during the phase transition, we have 

/ ~ 0.8 ( — ^Wt) dnuc, (1) 
VlOOMeV/ ' ^ ^ 
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where the factor Tc/lOOMeV scales the value of dnuc from a given critical temperature to 
T = lOOMeV, where / is defined. 

In the bag model we can approximate dnuc in homogeneous nucleation theory without 
curvature energy in terms of the bag constant, B, by 

/145MeV\2 fAB\ / a \3/2 
^"-"'H^V^J Ixj loOSSEVij (2) 

where L is the latent heat and a the surface tension. Lattice QCD calculations suggest 
that a < 0.24T3, Ref. [g finds a ~ 0.025r3 and L ~ OAB. With ~ 0.685^/^ this gives 
dnuc ^ 0.1m for > lOOMeV. 

We calculate the surface tension self-consistently within the multiple reflection expansion 
framework of the MIT bag model and get a ^ 0.0355^/"^ and L ~ AB, corresponding to 
dnuc ~ 0.1m. But we also self-consistently take into account the important curvature term 
and find a crude approximate formula for dnuc, 

, „ /145MeV\^ /45\ 
dnuc ^ IS ) (^jcm. (3) 

Note, that in this approximation the mean nucleation distance depends only on the bag 
constant (or T^). Thus the overall effect of including a curvature term is roughly a factor 
of 7 reduction of dnuc- With these estimates of dnuc in mind it hardly seems possible to get 
any significant deviations from standard nucleosynthesis, if homogeneous nucleation theory 
applies. 

However, homogeneous nucleation may not be the prime mechanism responsible for 
hadron formation. It is well-known from every day life that first order transitions are nor- 
mally facilitated by impurities, such as dust or ions and container walls. Think e.g. of a 
charged particle entering a cloud chamber, or boiling water where (half) bubbles are formed 
at the bottom of the pot. Such nucleations are denoted heterogeneous. Thus it is obvious 
to consider the effect of such impurities on the quark-hadron phase transition. Candidates 
are the topological impurities, like primordial black holes, magnetic monopoles, and cosmic 
strings or relic fluctuations from the electroweak transition. 

The idea that impurities could play a role in the cosmological quark-hadron transition is 
not new (see e.g. for a general discussion of the role of impurities in cosmological phase 
transitions), but normally it has been argued |^ that the presence of impurities would just 
act as extra seeds for nucleation sites together with the homogeneous nucleation sites and 
thus reducing dnuc- 

As demonstrated below, this picture is not necessarily correct — there are circumstances 
where impurities lead to a dramatic increase in dnuc-, maybe even to a value of several 
meters, resurrecting inhomogeneous Big Bang nucleosynthesis. We find that three limits 
with, respectively, small, medium and huge amounts of impurities, compared to the number 
of homogeneous nucleation sites, have to be considered. Certainly, if only very few impurities 
are present in the early Universe dnuc will be determined solely by homogeneous nucleation. 
In the other extreme, for a huge number of impurities, they will determine a dnuc smaller 
than the one obtained from homogeneous nucleation. So neither of these cases give large 
mean nucleation distances. 

Any impurity tends to lower the cost in energy when creating the surface of the new 
phase. Therefore bubbles of the new phase will nucleate easier around an impurity at a 
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smaller amount of supercooling. This means that shock spheres from nucleations around 
impurities start reheating the Universe earlier. And for some impurity number densities the 
whole Universe has already been reheated before the amount of supercooling has increased 
to a value where homogeneous nucleation sets in. In this case dnuc depends on the impurity 
number density, n, on the time when heterogeneous nucleation starts, and on the time when 
homogeneous nucleation becomes important. For a large range of parameters, dnuc can be 
very big. 

We will make this qualitative discussion much more quantitative in Section III, where we 
also present a self-consistent bag model calculation of homogeneous nucleation, including for 
the first time the important curvature terms. In Sec. II we describe the thermodynamics of 
the phase transition and the multiple reflection expansion of the MIT bag model, which we 
use to calculate the thermodynamical quantities. Section IV contains a general discussion 
and a summary of our results. 



II. THERMODYNAMICS AND THE BAG MODEL 



All the necessary thermodynamical information about the two phases is contained in the 
thermodynamical potential fi(T, fi, V, S,C). fi is the chemical potential, V the volume, and 
finite size corrections enter via the surface area, S, and the extrinsic curvature, C. 

n = -T\n Z, (4) 

where Z is the grand partition function. All other thermodynamical quantities can be 
derived from Q. For instance, the bulk pressure P is 

V^*^ / T,ii,S,C 

and the surface tension a 

(—] (6) 

Because of the huge photon to baryon ratio, we can neglect all chemical potentials and 
simplify the equations for Q when we want to calculate dnuc, since (fJ^/T) ~ 10~^ in the 
early Universe (see e.g. [|lOl)- (This is not true if one is interested in following the further 
evolution of baryon density contrasts.) 



A. The quark phase 



The quark phase consists of massless gluons, massless (or with masses of 5-10 MeV) u 
and d quarks, and a massive s quark, where we use a broad interval, rris/MeV G [50, 300]. In 
most calculations we have integrated numerically the s-quark distribution function; however 
for simplicity we neglect the small s-quark contribution to the cosmic density when calculat- 
ing the cosmological time-temperature relation from the Friedmann equations. The gluons 
contribute with 16 degrees of freedom, while each of the quarks has 12 degrees of freedom 
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(including antiparticles). In thermal equilibrium with both phases during the transition are 
(statistical weight in parenthesis) electrons(4), muons(4), neutrinos(6) and photons(2). The 
contributions from these particles cancel when we study pressure differences. 

The behavior of quarks is mainly governed by the strong interaction. It can, at least in 
principle, be described in the QCD model. But since QCD is a non-perturbative theory at 
low energy and therefore very difficult to handle, phenomenological models, like the MIT 
bag model, or lattice QCD calculations are used in practice. The basic idea of the bag model 
is that the quarks can be treated as an (almost) ideal gas confined to a finite region of space, 
a bag, by what we may think of as an external pressure, B, called the bag constant. In 
the multiple reflection expansion framework [jll| a statistical approach is taken and finite 
size effects of the bag are explicitly incorporated in the density of states formula Eq. (|l^) in 
terms of geometrical factors. 

What value to take for the phenomenological bag constant is at present very uncertain; 
it could even be temperature dependent. A lower limit on the bag constant of (145MeV)^ is 
set by the stability of nuclei, especially ^^Fe, relative to M(i-quark matter. There is no upper 
limit available from physical considerations. Lattice QCD calculations typically favor B^/^ 
above 200MeV. We use here a broad interval for the bag constant, i^^/'^/MeV G [145; 245]. 
It was shown, in the T = limit, by Farhi and Jaffe that a QCD coupling constant 



different from zero can largely be absorbed in a reduction of the bag constant. We assume 
this is valid at finite temperatures too and put = 0. 

The multiple reflection expansion was originally derived for a quark droplet embedded in 
a hadron gas. But Mardor and Svetitsky [|l^] have shown that it can also be used to describe 
a hadron bubble in a quark plasma, provided the signs of the volume and curvature terms 
are changed. The formulae below are for a quark droplet. The thermodynamical potential 
of the quark phase is then given by 

n, = Y.n, + Bv, (7) 

i 

where VLi is the thermodynamical potential for a single fermion (quark) or boson (gluon), 
given as 

Vti = -T\YiZi. (8) 
Here Zi is the corresponding grand partition function 

In = ± dk p{k) In (^1 ± e- (9) 

(upper sign for fermions, lower for bosons) where the density of states is obtained from the 
multiple reflection expansion 



+ fc(-]{(^ + ^]dS+...l (10) 



m j Js \Ri R 



g is the statistical weight of the particle. The volume term is universal and we recognize 
the finite size effects from the surface area, S, in the second term, and from the extrinsic 
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curvature, C, in the third term. Higher order terms are neglected. Assuming that the quark 
droplets are spherical, then V = AttR^/S, S = AttR^ and C = SirR. The coefficient functions 



fj{k/m) are for quarks 
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(12) 



/c has not been calculated from the multiple reflection expansion for massive quarks, but 
Eq. (|12D gives a perfect fit to shell model calculations for strangelets [|T5l. For the massless 
gluons we have |]T3| 



fs 



0, 



f 



c 



(13) 



Note, that these results for gluons have been conffimed by pure glue QCD lattice calculations 



B. The hadron phase 



We treat the hadron phase as an ideal pion gas, where we calculate the pressure from 
VLi using Eq. for bosons and only including the volume term in the density of states 
formula Eq. (0). This approximation is justified because the pions, with masses around 
the critical temperature 100-200MeV, are the lightest hadrons and semi-relativistic while 
all other hadrons are nonrelativistic and therefore contribute negligibly to the pressure for 
T < 200MeV. The pion pressure is much smaller than the bag pressure so even the pressure 
from fully relativistic pions will only increase the critical temperature one percent. A perhaps 
more serious problem is that pions are not "point-like" particles and excluded volume effects 
should be taken into account. The pions could in principle contribute to the surface term 
(and curvature term as well), we do not know. However, the pions have only 3 degrees of 
freedom, while each of the quarks have 12, so we would expect only minor corrections. 



C. Hadron bubbles in quark-gluon plasma 

The change in the thermodynamical potential, when a quark matter droplet of radius 
Rq is formed in a hadron gas, is 

^^q-droplet = PhYq + ^q{Rq)- (14) 

Because the photon-lepton gas is in thermal equilibrium with both phases, it does not 
contribute to Ai7. Only in pressure ratios is it of importance to include the photon-lepton 
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contribution. In the cosmological quark-hadron phase transition hadron bubbles are formed 
from a quark-gluon plasma. Mardor and Svetitsky |13] have shown that this situation can be 
described simply by inverting the system, i.e let Rg —Rh- This inversion affects only the 
signs of the volume and curvature terms and not the absolute value of any of the coefficients. 
Afl can then be written as 



bubble 



-P^Vh + nJ-Rh). 



(15) 



or 



An 



APR^ + AiraR^ - Sn-fR, 



(16) 



where we have suppressed the index h on the hadron bubble radius. If we relate the tem- 
perature and bag constant as 



then from the multiple reflection expansion we have 

37 

AP = Ph-Pg = P^ + B-Ps- 



(17) 



(18) 



where the pion pressure is (we treat the three pion types alike and use a mean mass of 
= 138MeV) 



32; 



27r2 

and the s quark pressure 

6x 



B duu^\n(l-e~^^ 



+0^ /x 



(19) 



(20) 



The critical temperature is found from AP = (the corresponding Xc lies between 0.67 and 
0.70). The surface tension, where only the massive s quarks contribute, is given by 



3x 



2n 



And finally the curvature coefficient 



a = (To = — B^^^ I du u i 1 arctan — 
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(21) 
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In (l + e-^^^l^ 



(23) 



7 



7 is always larger than zero because of the large positive contribution from the gluons, even 
though 7s can be negative depending on the mass of the s quark. So the curvature energy 
is always negative when making a hadron bubble. 

The latent heat per volume, released in a first order transition, is defined as 

(24) 

Eq. ( p!8D can easily be differentiated with respect to temperature and the latent heat calcu- 
lated numerically for fixed bag constant and strange quark mass. In Fig. |l] we have plotted 
the latent heat as a function of the strange quark mass. And we show the effect of including 
higher mesonic states up to 1 GeV (the higher mesonic states have much less influence on 
the critical temperature, because the bag constant contributes with about 95% of the pres- 
sure needed to maintain pressure equilibrium between the two phases. That is why we can 
safely disregard heavier mesons for T^. < 200MeV). It is evident that these states only have 
a minor influence on the latent heat. A good approximation is to use 

L ~ 45 (25) 

independent of m^. We adopt this approximation, but notice that some lattice QCD cal- 
culations suggest 1^ that the latent heat is much smaller, L ~ 0.45. We will discuss the 
possible discrepancy later. 



L = -T,. 



dAP 
dT 



III. NUCLEATION THEORY 

In classical nucleation theory the nucleation rate (number of nucleations per volume per 
time) is given by |I7| , p!8| 

p(T)=C(T)expf-^V (26) 



T 

where AF^ is the minimum work needed to create the smallest possible growing bubble, equal 
to the change in Helmholtz free energy, and exp(— AF^/T) is the probability of making 
a growing bubble. The pre-exponential factor depends in detail on the model by which 
the bubble growth is described. Csernai and Kapusta have calculated C(T) in the 
cosmological case (relativistic particles and almost zero baryon number), but without a 
curvature term, and get 

^(^) = 35/2^ ^4^22^3/2 (27) 

where Rc is the radius of a critical size bubble and ~ 0.7fm is the surface thickness of the 
bubble, has to be much smaller than R^, which it is since Rc ~ 50fm, typically. Finally 
the shear viscosity A is approximated by 

A ^ ^^T^ (28) 

where as is the strong coupling constant. Fortunately, as we will show, a temperature 
dependent pre-exponential factor like Eq. (pTl) is not important when we want to determine 
dnuc- Thus, the often used dimensional estimate C(T) ^ is sufficient for this purpose. 
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A. Homogeneous Nucleation 



This subsection is for a major part a straight forward generahzation of the work done by 
Fuller, Mathews and Alcock ||^ to take a varying pre-exponential factor and, more important, 
a curvature term into account. Their results (with the correction mentioned in pO|) are 
recovered in the limit of constant pre-exponential factor and vanishing curvature. 

The change in Helmholtz free energy is equal to the change in the thermodynamical 
potential in the limit of vanishing chemical potentials. Thus, 

AF{R, T) = APR^ + A'KaR^ - S7:-fR. (29) 

3 

The radius of a critical size bubble is found by putting dAF/dR equal to zero. From this 
equation we get two critical radii 



The smaller radius corresponds to a local minimum in the free energy (see Fig. 0), the 



larger to a local maximum. As pointed out in |T3| this local minimum has some serious 
consequences, if it is real and not just a shortcoming of the multiple reflection expansion of 
the bag model, since it means that small hadron bubbles show up even above the critical 
temperature. It has been argued that these hadron bubbles are unstable I^T], but at present 
it is unclear how to handle them. We will for now neglect this complication, keeping it in 
mind for later discussion. 

A linear expansion of AP about the critical temperature gives 

AP = Lr], (31) 

where L is the latent heat per volume, and the supercooling parameter for T <Tc is 

T —T 

V - (32) 

This expansion is valid for the small supercooling 77 <C 1 we deal with in the cosmological 
quark-hadron transition. We then write AF^ = F{R^,T) — F{0,T) as 

Certainly AFc can not be allowed to become negative if we insist on ignoring the minimum. 
This gives an upper limit on the supercooling parameter. 



Combining equations and ( P^D gives the nucleation rate as a function of the small 
supercooling parameter, valid around the critical temperature: 
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p{ri) = Cirj) exp 



Stt 



3 L^Tc?]^ 



(l - ?,hr] + (1 - 2hr]) 



3/2 



where 



(35) 



(36) 



We see that p{ri) equals zero as we pass the critical temperature from above and increases 
very rapidly with t]. As a result nearly all nucleation takes place at the lowest temperature, 
Tf, achieved during the supercooling phase. This means that we may generally approximate 
the pre-exponential factor with its value at Tf provided it is not too temperature dependent. 

It turns out to be fruitful to make a linear expansion of In p{ri) about the time tf corre- 
sponding to the temperature Tf 



\iap(t) ~ \np{rif) 



dlnp 




dr] 




dT 




dr] 


Vf 


dT 




dt 


if 



(t-tf), 



or 



(37) 



(38) 



p{t) ^ p{r]f) exp[-a{t - tf)]. 

This is a good approximation because of the steepness of the nucleation rate (cf. Eq. (|35D). 
We have implicitly assumed that the surface tension and curvature coefficient are indepen- 
dent of temperature. This is not true in the bag model. But there is nothing special at all 
about the critical temperature from the point of view of a and 7. And since the amount of 
supercooling is very small we treat them as constants and use their values at Tc. 

The temperature as a function of time is found by solving the Friedmann equations for 
a flat Universe consisting of a quark-gluon plasma, photons and leptons. It is, in the first 
order approximation we use, independent of the bag constant. 



T 



45 



1/4 



(39) 



where G is the gravitational constant, and the statistical weight Qq = 51.25 (for simplicity 
we here neglect the small contribution from the semirelativistic s-quarks). a is then 



a 



C'iVf) 



+ 



Civ) 

2 - 3br]f + (2 - brif){l - 2brifY/^ 



(40) 



where C = dC/drj. 

The mean nucleation distance is found from the total number density of nucleation sites 



f{t)p{t) dt, 

where f{t) is the unaffected fraction of the Universe. 



(41) 
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After a hadron bubble is nucleated a spherical shock front expands into the quark phase at 
a speed slightly above the sound speed, Vg — 3~^/^, and reheats it to the critical temperature, 
which prevents any further nucleation in a region crossed by one or more shock fronts [Q. 
We use the formula given by Guth and Tye p2| for the fraction of space not yet passed by 
one or more shock fronts (a slightly different approach can be found in [0), 



f{t) = exp 



- / dt'f{t')p{T{t'))V{t',t) 



(42) 



where tc is the time corresponding to temperature and V{t',t) in this context is the 
reheated volume at time t caused by a bubble nucleated at time t'. V(t',t) can be written 

as 

Vit',t)^fv!it-tr. (43) 

For the short timescales involved here, we can safely neglect the cosmic expansion. This 
recursive formula for f{t) can be simplified to a step function by a linear expansion, giving 

^ I J - iL dt'pinmt - t'f t < t^f _ ^^^^ 

The time tf corresponds to the maximum supercooling rjf, where we demand the whole 
Universe to be reheated. The approximation succeeds because of the very rapid rise of the 
nucleation rate. If we introduce the approximation for the nucleation rate, Eq. (|38|), the 
integral in Eq. (^if) can be analytically evaluated and 



SttV 



^ ^ \ t>tf ^ ^ 

where terms exp [—a{tf — tc)] have been neglected compared to terms exp [—a{tf — t)], be- 
cause of the steepness of p{t). Demanding the whole Universe to have been reheated at tf 
gives an equation for the maximum amount of supercooling achieved during the transition 

-^PiVf) = 1. (46) 
a 

It has an approximate solution given by 

0.4 (t3/2 



^3/2 



^0-48-rT7^\/l-3&r//. (47) 

A finite 7 approximately halves the critical amount of supercooling, since it turns out that 
hr\f ~ \. The behavior of r^/ is shown in Fig. |^. Notice, that 77/ scales with 1/L. 

The approximate solution in Eq. ( ^7] ) ignored the details of the pre-exponential factor, 
C . To see that this is justified, neglect for a moment the curvature term. It is then easy to 
show that 
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-1/2 



with all dimensional quantities in MeV-units. For the C{ri) in Eq. (|27| ) C'/C = —l/rjf. 
Because 

3L^r]j rjf 

when rjf is approximated by Eq. (^Tf), we can safely ignore the C'/C term. Thus a pre- 
exponential factor of the form Eq. (^TD does not significantly influence the maximum amount 
of supercooling neither directly, not even changes in C of several orders of magnitude, nor 
through its derivative. This result is independent of bag model parameters. The same qual- 
itatively picture holds when the curvature term is included, if we do nothing but substitute 
Eq. (^) for i?c. Notice however, that this may not be valid since the pre-exponential factor 
is derived without taking a curvature term into account. 

The mean nucleation distance is found by substituting Eqs. (^) and ( pH]) into Eq. (^Tj) 
and integrating, 

d„,,= (167r)i/3!^, (50) 
a 

where terms like exp [— a(t/ — tc)] have been neglected. It should be emphasized that this 
is also a general result independent of bag model parameters. Only a is model dependent. 
With a from Eq. (EOf), neglecting the C'/C term, dnuc becomes 



dr,,,, ~ 1.3 X 10^ 



-MeV^ 



cr 



^/ (51) 



-I 

cm. 



Eq. d^) is found by inserting the first approximation for rjf from Eq. (^) with 7 = 0. For 
7 7^ one finds in the same manner from the second approximation in Eq. (^) 

- 1.8 i^^)' if) i^fMBU^)' (52) 

but 7 is very insensitive to varying s quark mass because of the large contributions from the 
massless species. Thus the last term is always close to unity and we recover Eq. (|). Note, 
the mean nucleation distance is also inversely proportional to the latent heat. 

In Figure ^ we have plotted the mean nucleation distance as a function of the strange 
quark mass for a small and a large bag constant with and without the curvature term. 
Curves for intermediate bag constants lie between those two curves. This has been done 
by solving Eq. ( |iBD numerically and inserting in Eq. ( pUD for dnuc- For comparison the two 
approximative solutions, Eqs. (0) and (^ have been plotted, too. If we look at Eq. ( pOD we 
see that for increasing a, i.e. decreasing amount of supercooling, then the mean nucleation 
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distance decreases. This is because a larger a means a more sudden nucleation and more 
bubbles are needed in order to reheat the whole Universe at r]f. Inclusion of a curvature 
term lowers the energy barrier and thereby the maximum amount of supercooling, resulting 
in roughly a factor of 7 reduction of dnuc- 



B. Heterogeneous Nucleation 

In a first order phase transition the presence of impurities lowers the energy barrier and 
thereby reduces the maximum amount of supercooling achieved during the transition. We 
assume that a number density, n, of impurities were present in the Universe prior to the 
quark-hadron transition. Candidates could be primordial black holes, magnetic monopoles, 
cosmic strings, or relic fluctuations from the electroweak transition. Especially the latter 
candidate is interesting because the horizon distance at the time of the electroweak transition 
~ 0.3 cm [t ~ 10~^^s, T ~ lOOGeV) has expanded three orders of magnitude to ~ 3 meters 
at the time of the quark-hadron transition. Further we assume that nucleation around all 
impurities takes place at a time ti, where tc < ti < tf, and that the bubble expansion is 
similar to the homogeneous case. 

The combined nucleation rate can then be written as a sum of the nucleation rate around 
impurities and the homogeneous nucleation rate, 

Pcomit) = n5{t - t,) +p(r^^)e-"(*/-*), (53) 

where 5 is Dirac's delta function. The unaffected fraction of the Universe can still be eval- 
uated in the same approximation leading to Eq. (^) although the approximation becomes 
slightly poorer. 



1 _ ^JlA {n(t - Uf +-^p(,^^)e-"(*/-*)) t < tf 



^ ^ \0 t>tf ^ ^ 



and the equation for rjf now reads 



s 



The number density of nucleation sites is found from Eq. (PD, 



-3 , PiVf) 47rt;Xr/j)2 



a 



- 3/ + 6;^ - 6 



(56) 



where z = a(tf — ti). Terms like exp{—a(tf — t^)) have been neglected. We can write the 
time difference tf — ti as 



tf-t, = {l-i){tf-t,), (57) 



where < ^ < 1. 
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Eq. (p6D can be simplified by introducing Eq. ( p5D and keeping only the term in the 
parenthesis. Since a(tf — t^) ~ 600 and we demand > 10(32;^ — 6z + 6) then we must 
have ^ < 0.95. This is no serious limitation, as we will see. And we hereby decouple the 
heterogeneous and homogeneous contributions to the number density of nucleation sites. 
Eq. (E6D now reads 



dnuc -n + ^Tivl ~ n + — — . (58) 



~ n H 



In the limits of vanishing and huge n, we recover the expected results, dnuc = dnucAom (Eq. 
(0)) and dnuc = n~^^^, respectively. 

Because of the steepness of the homogeneous nucleation rate all nucleations take place 
within a 3% interval in r], just below rjf. This means that either homogeneous or heteroge- 
neous nucleation determines the mean nucleation distance. Only in a very narrow range are 
there significant contributions from both. 

In Figure |^ we have plotted dnuc versus n for fixed bag constant and different s quark 
masses, in the most optimistic case where ^ <^ 1 . Notice the step function behavior and the 
broad range, more than two orders of magnitude, of impurity number densities, resulting 
in mean nucleation distances above one meter. The largest mean nucleation distance for 
given parameters is found when — fj/^hom- This is a scenario where the shock waves from 
the nucleated hadron bubbles around the impurities have reheated the whole Universe just 
before homogeneous nucleation becomes important. The corresponding n can be found from 
Eq. (^5D by neglecting the homogeneous term, 

d„„,_ = ,r'/» . 6(1 - (3^) (^)' m. (59) 

The corresponding maximal value of /, comoving at T = lOOMeV, is found by combining 
Equations (^9]) and (|l|) into 

If ^ < 1/2 and with a typical value of r^j ~ 3 x 10~^ (for L ~ AB; otherwise rjf oc 1/L) one 
can have significant nucleosynthesis effects (/ of some meters) from impurities for a broad 
range of n. For ^ ^ 1 the approximations used in the derivations above break down, but in 
this regime the presence of impurities do not lead to large nucleation distances anyway. 

The surface tension rises steeply with increasing quark mass for small masses. For 
realistic u and d quark masses rriu = 5MeV and nid = lOMeV the surface tension 
increases with about 30% when the s quark contribution is near its maximum. And dnuc,max 
increases almost 70% to about 11 meters. If we artificially double the surface tension to 
mimic eventual hadron contributions, we gain more than a factor of 3 on the maximum 
mean nucleation distance. 



IV. CONCLUSION 

We have calculated the mean nucleation distance in a homogeneous nucleation scenario 
where a curvature term contributes to the surface effects. The framework of the multiple 
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reflection expansion of tlie MIT bag model has been used to calculate the surface tension 
and curvature term self-consistently. Inclusion of a curvature term reduces dnuc by about a 
factor of 7 to less than 2cm. 

Furthermore, we have shown that the presence of impurities, for a broad range of den- 
sities, can result in mean nucleation distances above one meter. Such large distances are 
needed if baryon diffusion should not have smeared out generated baryon inhomogeneities 
before the epoch of nucleosynthesis. It must be emphasized that this conclusion is qualita- 
tively independent of the model used to describe surface effects. We have only assumed that 
impurities are present, that they lower the energy barrier and that the reheating mechanism 
is the same regardless of how the hadron bubbles nucleated. Thus, similar effects could be 
important in other first order transitions as well (see e.g. Ref. 0). 

Unfortunately there are no reliable theoretical expectations for the nature or density of 
impurities, but we note again the coincidence, that the horizon distance at the electroweak 
phase transition has expanded to the interesting range of a few meters at the quark-hadron 
transition. 

Eq. ( |T7| ) for rjf in the homogeneous nucleation scenario tells us how uncertainties in the 
different parameters influence dnuc, which depends linearly on rjf. We believe the curvature 
coefficient is well-known, because the major contributions come from gluons and the two 
massless or very light quarks. A significant surface tension contribution from the pions 
would raise r]f and thus result in a larger dnuc,max and an even broader interval of impurity 
number densities, where the mean nucleation distances is above one meter. 

We have throughout used the self-consistently calculated latent heat of 45, even though 
some recent lattice calculations suggest that the latent heat is about ten times smaller. 
A temperature dependent bag constant that increases with increasing temperature may 
reproduce such a result, but no theoretical predictions of such a temperature dependence 
exist. However, rjf, dnuc, and / are all inversely proportional to the latent heat. Thus, a 
tenfold decrease in L results in a tenfold increase in Tjf, dnuc, and /. 



In Ref. [20 1 the spectrum of nucleation site separations is calculated. A duality between 



the nucleation sites and large baryon concentration sites is assumed and the effect on nu- 
cleosynthesis calculated. If impurities play a major role, the nucleation site spectrum would 
look different. For randomly distributed impurities the spectrum would be broad. A sort of 
regularity in the impurity distribution would result in a very narrow distribution around the 
mean nucleation distance. Such an approximation is often used in actual inhomogeneous 
nucleosynthesis calculations. 

A serious problem with the bag model in the form used in this paper is the local minimum 
in the free energy at about 10 fm, even above the critical temperature \T3\. One might think 
about removing it by introducing further corrections in the density of states related to zero 
point energy, color singlet restrictions, etcetera, but no recipe exists at the moment. We note 
however, that such effects will be most important for small radii, and since the maximum of 
the free energy barrier typically lies beyond a radius of 50 fm in the results presented above, 
the conclusions of the present investigation are not likely to be affected much. 
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FIGURES 



FIG. 1. The latent heat as a function of the s quark mass for 51/4 = l45MeV. The uppermost 
curve includes only pions, whereas the other curves from above include mesons up to masses of 
547MeV, 892MeV and 1020MeV, respectively. 

FIG. 2. The free energy as function of bubble radius for different temperatures in MeV at 
B^^^ = 145MeV and nis = 150MeV, with curvature. The critical temperature gives the solid 
curve. The dashed curve is for the maximum supercooling, rjf, where almost all hadron bub- 
bles are nucleated with a radius, i?+ ~ 55fm. The dotted curve is still for r]f but with 7 = 0. 
The dash-dotted curve corresponds to an even larger supercooling where the local minimum at 
i?_ ~ lOfm can no longer be ignored. 

FIG. 3. The maximum supercooling, rjf, versus m^. Solid curves are for = l45MeV, the 
upper for 7 = 0. Dashed curves are for 51/4 = 245MeV, the upper for 7 = 0. The dotted curves 
correspond to the approximation in Eq. (]47|). 

FIG. 4. The mean nucleation distance as a function of the s quark mass. The solid curves are 
for = 145MeV, the uppermost for 7 = 0. The dotted curves for B^/^ = 245MeV. Dashed 
curves are the corresponding approximative expressions in Eqs. (^) and (^). 

FIG. 5. The mean nucleation distance as a function of the impurity number density for different 
s quark masses at S^/^ = 145MeV. From left to right the masses are lOOMeV, ISOMeV and 
200MeV, respectively. 
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